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ABSTRACT 

Two  prediction  models  for  graphical  prognosis  are  developed  and 
evaluated.  The  primary  feature  of  these  models  which  represents  a  de- 
parture from  previous  graphical  methods  is  the  inclusion  of  the  ter- 
rain-induced vertical  velocities.  These  vertical  motions  may  in  turn 
produce  vorticity  changes,  and  thus  height  changes. 

The  first  model  presented  is  comprised  of  two  separate  parts;  the 
first  part  is  a  prediction  of  the  5>00-mb  surface  and  the  second  part  a 
prediction  of  the  1000-mb  to  500-mb  thickness.  The  two  parts  may  be  com- 
bined to  produce  a  surface  (1000-mb)  prognosis.  This  model  failed  to  pro- 
duce satisfactory  results,  and  hence  a  second  model  was  developed. 

The  second  model  presented  follows  somewhat  the  development  of  the 
first,  but  avoids  an  assumption  relating  the  surface  and  500-mb  winds 
that  is  believed  to  have  contributed  largely  to  the  original  failure. 
This  model  operates  on  the  current  1000-mb  and  1000-mb  to  5>00-mb  thick- 
ness charts;  and,  utilizing  a  fraction  of  the  5>00-mb  height  change  as 
predicted  by  any  other  established  method,  produces  a  surface  (1000-mb) 
prognosis.  The  model  was  tested  in  two  ways;  one  with  the  terrain  factor 
included  in  the  prediction  equation  and  one  omitting  the  terrain  factor. 
Both  methods  produced  satisfactory  prognoses.  However,  inclusion  of  the 
terrain  factor  consistently  produced  the  better  results. 

The  writer  wishes  to  express  his  appreciation  for  the  assistance 
and  guidance  of  Professor  George  J.  Haltiner  of  the  U.  S.  Naval  Postgrad- 
uate School  in  the  preparation  of  this  paper. 
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PREFACE 

Graphical  methods  of  numerical  prognosis  have  recently  received 
widespread  attention  since  the  development  by  Fjortoft  {jfj  of  a  predic- 
tion model  for  the  500-mb  surface.  Fjortoft  provided  a  graphical  tech- 
nique for  integration  of  the  barotropic  vorticity  equation.  Despite  the 
relative  simplicity  of  the  method  and  the  rather  restrictive  assumptions 
inherent  in  its  development,  tests  by  the  Air  Weather  Service  of  the 
United  States  Air  Force  [2]   indicate  accuracy  comparable  to  that  achieved 
by  a  trained  forecaster  using  conventional  techniques. 

Following  the  success  of  Fjortoft's  model,  several  other  more  in- 
volved graphical  models  have  recently  been  devised.  Although  none  as  yet 
have  been  subjected  to  the  degree  of  testing  that  Fjortoft 's  model  has, 
preliminary  results  are  encouraging. 

The  graphical  methods  devised  thus  far  appear  to  have  been  least 
successful  in  mountainous  terrain  areas.  The  reason  for  this  probably 
lies  in  the  simplified  boundary  condition,  CO  —  0  at  p0,  which  has  been 
assumed  in  the  previous  models. 

It  is  the  purpose  of  this  paper  to  attempt  to  incorporate  terrain 
effects  into  graphical  techniques  of  prognosis.  The  development  here  is 
somewhat  similar  to  the  treatments  of  Estoque  [S]   and  Reed  fiij  . 
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CHAPTER  I 
FIRST  PREDICTION  MODEL 
The  first  model  is  comprised  of  two  partsj  the  prognosis  of  the 
500-mb  surface  and  a  separate  prognosis  of  the  1000-mb  to  £00-mb  thick- 
ness. The  two  prognoses  may  be  combined  for  a  resultant  surface  ( 1000-mb) 
prognosis.  The  development  here  parallels  that  of  Estoque  J^3J» 

1,  Development  of  the  prediction  equation  for  the  5>00-mb  surface. 
The  "vertical  velocity"  is  assumed  to  vary  thusly: 


W(xj3(ftt)=-Kg/?(y5-VHX-f)4^^At) 


■? 


Vf.-v 


z"1 


(1) 


The  first  term  on  the  right  represents  the  contribution  due  to  the  oro- 
graphic effect  and  is  assumed  to  vary  linearly  with  pressure  (p/p  ), 
from  a  maximum  at  p  to  zero  at  the  top  of  the  atmosphere,  k  is  a  scalar 
multiplier  with  possible  values  ranging  between  zero  and  unity.  It  is 
assumed  to  equal  unity  for  this  development  and  hence  will  not  appear 
in  subsequent  equations  in  this  chapter.  The  second  term  on  the  right 
is  the  contribution  due  to  the  non-orographic  effects j  and  is  assumed 
to  vary  parabolically  from  zero  at  P=P0  to  a  maximum  of  0Om  at  p-p^. 
Inherent  in  this  distribution  is  the  assumption  that  p^  (to  be  taken  as 
^OOmbs)  is  a  level  of  non-divergence  where  there  are  no  orographic  effects. 
The  non-orographic  contribution  assumed  here  varies  slightly  from  Estoque' s 
assumption  of  a  sinusoidal  variation  of  vertical  velocity  j~3J  • 
Differentiating  (1)  with  respect  to  pressure  yields 

ap---ft-(y«-yH)-2e°»5£%f  .     (2) 


Evaluating  (2)  at  p^   gives 

The  following  form  of  the  vorticity  equation  is  assumed: 

This  form  of  the  vorticity  equation  neglects  the  vertical  advection  of 
vorticity,  the  vortex  tube  term,  and  J°7  ,  Since   /Sio  -  ~ V  ■  V  ,  Eq. 
(3)  shows  that  the  only  divergence  considered  at  ^OOmbs  ±s   due  to  ter- 
rain-induced vertical  velocities. 
Combining  (3)  and  (k)   yields 

In  order  to  obtain  conservation  equations  utilizing  only  £00 -mb 
data,  the  surface  wind  must  be  eliminated  from  Eq.  (5).  This  will  be 
done  by  assuming  that  JTg-aV,  where  a  is  a  scalar  multiplier.  Then 

(5)  may  be  written 

-+V-V1=-^V-VH    .  (6) 

Taking  mean  values  of  f  and  /°     for  the  area  to  be  prognosticated, 

(6)  may  be  written 

With  the  assumption  of  geostrophic  vorticity,   this  becomes 

|^W-vffV£++  H-a^LH)=0      ;or 


^+y-v)(fy2+f+^-H)=o 


dt 


(8) 


as  f  and  agf/£H/p0  are  invariant  with  time. 

Introducing  the  finite  difference  formula  for  vorticity  into  Eq. 
(8)  gives 

r4 


(* 


s+y-v)(^(^)^  +  ^LH 


o 


(9) 


Multiplying  through  by  factor  fd2/W^g>  which  is  invariant  with  time, 
and  expanding,  produces 


at 


2  »* 


(a) 


V 


u 


C8^  +  £^+fifdJfjSH    + 

(b)  (c)  (d) 

2^M+V.^7f  +  V.^>««£SH  =  0    (10) 


Considering  terms  of  Eq.  (10)  separately: 

(a)  This  term  may  be  written  9/9t(Z  -  Z),  as  other  terms  are 
invariant  with  time. 

(b)  The  factor  ljm2g/fd2  may  be  taken  outside  of  V  operator,  as 
variation  is  small  compared  with  the  variation  of  (Z  -  Z) . 

(c)  Following  Fjortoft  \l\,   this  term  may  be  written  as 


\AV«J(9)>  where 


J(<f)  = 


'  S£  d  smf  gosh> 


A 


m*3 


df 


(d)  Taking  mean  value  of  f  for  area  to  be  prognosticated,  this 
term  may  be  written 


—  •3 — 


v-v 


a 


H  wzf> 


H 


Rewriting  (10)  with  above  changes  gives 


Hmfcf0 


-O    .  (11) 


With  the  notation 

B-4^2H   > 

Eq.  (11)  becomes 

81:       —  *"  J 

V— 

As  B  and  J(cp)  are  invariant  with  time,  (12)  may  be  written 


(12) 


or,  in  Jacobian  form, 

^  (Z-Z  +  Jcp)  +  B)  +tJ(S,  Z~Z  +Jw+B)  =  °  ^     (1W 

which  is  the  desired  prediction  equation.  However,  in  using  the  above 
equation,  the  advection  would  be  accomplished  with  the  current  500-mb 
wind,  which  obviously  may  vary  considerably  with  time,  especially  if  a 
prognostic  interval  of  12  or  2k   hours  is  employed.  In  order  to  arrive 
at  an  advective  field  which  is  more  conservative  with  time,  the  Jacob- 
ian identity  j(X,Y)=£j(X-f-I,  Y)  may  be  employed.  Applying  this  to  (Hj.) 
gives 

Thus  the  advection  is  accomplished  with  a  much  more  conservative  field, 
("Z+  J(f)  +  B),  Z*  being  the  only  quantity  that  varies  with  time,  and  this 
much  more  slowly  than  Z. 

As  a  preliminary  test  of  results  of  the  model  developed  herein, 
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Fig.  1  was  prepared.  An  initial  west-east  zonal  flow  of  $0   knots  was 
assumed  at  ^OOmbs.  The  surface  wind  was  assumed  to  be  equal  in  direc- 
tion and  J  the  magnitude  of  the  500-mb  wind  (i.e.;  a  =  ^) .  Eq.  (l£) 
was  used  and  a  prognostic  chart  for  a  time  of  six  hours  was  prepared, 
using  height  changes  as  approximated  by  the  first  term  of  Fjortoft's 
series  [l].  The  resultant  prognostic  chart  shows  a  trough  in  the  lee 
of  the  Rocky  Mountains,  as  would  be  expected  under  the  circumstances. 

2.  Testing  of  the  prediction  equation  for  the  ^OO-mb  surface. 

The  first  step  in  the  employment  of  Eq.  (15)  was  the  construction 
of  the  Z*  chart.  This  chart  was  prepared  by  conventional  graphical  methods, 
using  a  grid  distance  d  of  HKkm  as  suggested  by  Estoque  [3j  •  Next,  the 
contours  of  J((p)  corresponding  to  a  grid  distance  of  ICKkm  were  added 
graphically  to  the  Z  chart.  Assuming  a  value  of  ^  for  a,  the  quantity  B 
was  evaluated  as  1.5  x  10  ,  assuming  standard  values  of  the  meteorolog- 
ical parameters,  neglecting  variations  of  m  ,  and  taking  d  as  10-^km,  as 
above.  (Variations  of  m  are  small,  as  a  conformal  projection  with  stand- 
ard latitudes  of  30N  and  60N  was  used,  and  prognostication  was  confined 
to  middle  latitudes).  A  terrain  contour  chart  of  North  America  was  multi- 
plied by  the  above  constant  and  isolines  of  the  quantity  B  were  drawn. 
The  B  chart  was  then  added  to  the  Z  +  J(f)  chart  to  produce  the  advective 
field  (Z*  -f  J(f )  -|-  B)  .  The  Z  chart  was  then  subtracted  from  the  advective 
field  to  produce  the  (Z  -  Z-f  J(^ )  +•  B)  field,  or  the  conserved  quantity 
to  be  advected. 

The  above  two  charts  were  superimposed  in  register,  and  isolines  of 
(Z  -  Z  +  J(<f)  +  B)  advected  geostrophically  with  the  (Z  +  J(<f)  +  B)  chart 
for  the  chosen  prognostic  interval  of  12  hours. 


Fig.  1 


The  prognostic  12 -hr  height  changes  may  have  been  recovered  by  use 
of  Fjortoft's  series  [l]  , 

AE=A(z-Z+Jc<f>+B)  +  2(H-Z+0«f>+B) 

and  compared  with  observed  values  of  height  changes,  but  for  this  test  a 
more  direct  method  of  evaluation  was  used.  The  quantity  ("Z"  -  Z  +  J((p)  +  B) 
was  advected  for  the  12 -hr  period  and  compared  with  the  observed  values 
of  the  same  quantity  as  computed  from  the  observed  chart  at  the  end  of  the 
prognostic  period. 

The  results  of  four  tests  using  the  above  techniques  and  assumptions 
were  generally  poor.  Indications  were  that  in  some  isolated  regions  the 
inclusion  of  the  terrain  factor  gave  better  results  than  the  use  of  the 
prediction  equation 

^  (z -Z  +  J(T))  +  £-  J  (H  +  J«f),  2  -Z  +J«f>) -  O.  (16) 

(Except  for  the  different  vertical  velocity  distribution  assumed,  this 
equation  is  essentially  the  same  as  Estoque's  model  [YJ ) .  However,  these 
regions  of  improved  results  were  not  constant  in  space  nor  were  they  pre- 
dictable, and  the  areas  of  poorer  results  generally  equalled  or  exceeded 
the  improved  regions,  resulting  in  generally  poorer  prognoses  by  inclusion 
of  the  terrain  factor  than  by  ignoring  it. 

It  was  reasoned  that  the  choice  of  a=  ^,  i.e.,  that  the  surface 
wind  is  equal  in  direction  and  \   the  magnitude  of  the  500-nrt)  wind  (in 
connection  with  the  terrain  factor),  may  have  been  largely  responsible 
for  the  failure  of  the  method.  Therefore,  another  technique  was  develop- 
ed using  a  equal  to  either  \3   zero,  or  -§-,  as  determined  in  the  follow- 
ing manner.  The  current  5>00-mb  and  850-mb  charts  were  superimposed,  and 
areas  delineated  where  the  £00-mb  contours  were  approximately  the  same 


direction  as,  at  right  angles  to,  or  opposite  in  direction  to  the  Q$0-mb 
contours.  The  850-mb  chart  was  selected  as  the  chart  most  nearly  represent- 
ing the  surface  in  the  majority  of  the  mountainous  areas.  Then,  in  pro> 
ducing  the  (Z  -f  J(y> )  +  B)  chart,  the  contours  of  B  as  drawn  for  the  first 
technique  were  added,  considered  zero,  or  subtracted  in  the  areas  indicat- 
ed. The  subsequent  steps  in  this  technique  were  as  before. 

Results  of  using  the  new  technique  were  slightly  better  than  for  the 
assumption  that  a  is  everywhere  ^.  However,  improvements  realized  were  in- 
sufficient to  justify  use  of  the  method.  A  prime  difficulty  in  applica- 
tion  was  the  determination  of  appropriate  values  of  a,  which  appeared  to 
vary  rapidly  over  small  areas,  especially  in  situations  where  closed  lows 
were  present  at  the  5>00-mb  or  850-mb  levels.  An  extensive  search  was  made 
to  locate  a  synoptic  situation  where  the  differences  between  500-rnb  an<x 
85>0-mb  wind  fields  were  clear-cut,  in  an  attempt  to  evaluate  the  technique 
objectively.  Such  a  situation  was  located,  and  results  of  the  testing  in- 
dicated that  whereas  the  sign  of  a  was  correctly  chosen  by  the  above 
technique,  the  magnitude  of  a  appeared  to  vary  considerably  from  the  as- 
sumed values  of  \t   zero,  and  -g.  As  a  consequence,  the  prognostication 
was  not  anpreciably  better  than  that  using  Eq.  (16),  and  indications  were 
that  even  if  the  sign  of  a  were  predictable,  the  magnitude  would  vary  to 
such  a  degree  that  prognostications  ignoring  the  terrain  factor  would  be 
more  reliable. 

3.  Conclusions  as  to  prognosis  of  the  |?00-mb  surface. 

No  improvement  over  Estoque's  prediction  model  Pfj  is  realized  by 
inclusion  of  terrain  effects  on  vertical  motion  by  methods  developed  in 
this  paper.  It  is  believed  that  the  major  reason  for  the  failure  lies  in 
the  difficulty  of  expressing  the  surface  wind  in  terms  of  the  500-mb  wind 
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in  the  500-mb  conservation  equation. 

km   Development  of  the  prediction  equation  for  the  1000-  to  5>00-mb  thick- 
ness, 

A  detailed  development  of  the  prediction  equation  will  not  be  pre- 
sented here.  The  development  follows  closely  the  methods  of  Section  1  of 
this  chapter,  so  a  listing  of  the  major  assumptions  and  the  final  equa- 
tion will  be  shown. 

Assumptions: 

(a)  The  vertical  velocity  distribution  is  as  described  in  Eq.  (1), 

(b)  Temperature  changes  are  adiabatic.  The  following  form  of  the 
adiabatic  thermal  equation  is  employed: 

at 

(c)  Vr  'V^f  and  Vn'V-f  are  relatively  small  and  are  neglected. 

(d)  Similar  assumptions  as  those  used  in  progressing  from  Eq.  (10) 
to  Eq.  (11)  in  Section  1  of  this  chapter. 

With  the  notation 


£.(££) +v-v"s£  +GO<T=0  (17) 


r         qaJd/gH       ad     K   _ ,   ,  Sifd 


Z 


the  prediction  equation  evolves  as 


(18) 


or,  in  Jacobian  form; 


Due  to  the  failure  of  the  500-mb  model  to  produce  an  acceptable  prog- 
nosis, no  attempt  was  made  to  test  Eq.  (19).  Rather,  a  different  develop- 


ment  was  attempted  in  which  the  assumption  J/g  =  aV  was  not  necessary,  as 
it  is  believed  that  this  assumption  was  the  major  cause  of  failure.  The 
new  development  is  presented  in  Chap.  JI» 
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CHAPTER  ZT 
SECOND  PREDICTION  MODEL 
This  model  is  designed  to  be  applied  in  conjunction  with  a  separate 
3>00-mb  prognosis,  such  prognosis  to  be  made  by  any  method  the  user's  ex- 
perience deems  the  most  accurate.  Fjortoft's  method  [YJ  is  suggested*  The 
development  herein  parallels  that  of  Reed  \kj . 

1.  Development  of  the  prediction  equation. 

The  vertical  velocity  distribution  of  Eq.  (1)  is  also  assumed  for 
this  model.  Differentiating  (1)  with  respect  to  pressure  yields 

,        (20) 


V 

US* 

vh) 

-Zoj^ 

-p5f 

and 

evaluation  at  lOOOmbs 

i  gives 

~  — 

K9^ 

fo 

vh)- 

(21) 

Assuming  that  temperature  changes  are  adiabatic,  the  following  form 
of  the  adiabatic  thermal  equation  is  employed: 

By  integration  from  p  to  p^,  (22)  becomes 

V 

|£+VS- Vh  =  crfcocip    .  (23) 

The  surface  wind  may  be  used  for  the  advection  in  lieu  of  a  suitable 
mean  wind  for  the  layer  by  considering  the  wind  hodograph  to  be  linear 
and  the  thermal  wind  relationship  to  hold,  i.e.  : 
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Evaluating  the  right  side  of  (23)  by  use  of  (l)  leads  to 

crjf  licLp  =  -  cr  Kg  a  (Vs  •  VH)  f  j  \ 


-f  crco 


r, 


z 


-JK£*)>- 


Performing  indicated  integration,  and  noting  that  P0=  2p^,  gives 


cr     codf  --crKg^(\/$-VH) 


y^fV 


£ 


3(f0-f5) 


which  finally  reduces  to 

cr]  (Ldp  =  -lo-Kg/g(y,-7H)fs+|.<rncoB1  . 

Substituting  from  Eq.    (21;)   into  Eq.    (23)  produces 


^ 


(2U> 


3h  3 


it+Zs^h+fcrK3^(vs-7H)i|=4crfscuni  <*> 


Evaluating  Eq.   (U)   at  p0  yields 


at 


'  'o 


(26) 
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(27) 


assuming  that  Vg  =  V0#  This  assumption  is  used  throughout  this  develop- 
ment. Combining  (26)  and  (21)  produces 

Multiplying  (2£)  by  factor  3f /a~  P^  gives 

which,  upon  substituing  in  (27).  yields 


(23) 


3t 


U, 


'  ~-~     «v-v». 


3+ 


cr 


1°; 


ah 
3 


f+y''yh+|o-K^(vs.7Hjfs 


(29) 


J 


Collecting  advective  terms  and  time  derivatives,  and  again  noting  that 
pQ—  2p^,  (29)  may  be  written 


M**+&h)---*-zk+&h+i+^H 


Making  the  geostrophic  approximation  for  J0  leads  to 

3-f 


(30) 


5.  [  H  m-g  /- 

^2  l*o  W +crfr2 


h 


vs-v 


id2 


3f 


(?o-Zo)+^h+4  + 


(Tf, 


K3-f/lH' 


which,   after  multiplying  by  factor  fd  /Iiflrg,  yields 


5 
2/i 


21? 


(3D 
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9. 

at 


Z0-?o)  + 


3+V 


Mmzgcrf£ 


-V,-\7 


2  —  *  Z 


0/ 


Hm3crp5 


^fo 


(32) 


where 


Y 


J< 


(f) 


o 


m2g 


following  Fjortoft  £lj .  With  the  notation 


r_  IIK-fV/% 

8  m*  ft 

(32)  may  be  written 

3 

at 

?0-?0+Ch  +J(^ 

and 


C  = 


d<p 


as  J ((f)   and  F  are  invariant  with  tine.  This  is  the  prediction  equation, 
which  may  be  written  in  Jacobian  form  as 


ax 


_.    2 


Ho-Zo  +  OhJc^+F  =-f  J  Zo,Z0-2otCh+Jc^)  +  F 


(33) 


However,  if  (33)  were  to  be  used  in  prognosis,  the  advection  would 
be  accomplished  with  the  current  1000-mb  chart,  which  may  vary  consider- 
ably over  a  prognostic  period  of  12  or  21;  hours.  In  order  to  arrive  at  a 
more  conservative  advective  field,  the  Jacobian  identity 

J(xJY)  =  J(x  +  Yjr) 

may  be  employed.  Applying  this  to  (33)  leads  to 

ar 


-  ^-j[Zo+Ch+J^  +  F,Zo~Z0-fCh-fJc^4-F]  ? 


(3U) 


lii 


which  accomplishes  the  advection  with  a  more  stable  field  than  Z  . 

To  recover  height  changes,  the  local  change  in  the  advected  quantity 
over  the  prognostic  period  may  be  written 


A 


Z0-Z0  +  Ch  ■+ J(q>)  +  F 


=  -A 


or 


Z\[zo-?D+Chj  =  -A 


(35) 


as  F  and  J  (if )  are  invariant  with  time.  With  the  relationship  h  =  Z?  -  ZQt 
(35)  may  be  written 


A 


?o-H0+Cfs-£20"j  =  -A 


Solving  for    A  Z0  yields 


A20  = 

Utilizing  Fjortoft's  method  [l]  to  evaluate  A^0  gives 


A^0  - 


C  + 


A  -V   CAZ5 


Substituting  (38)  into  (37)  yields 

1   r 


A2A  = 


c-t- 


A  +  -£■  A 

C-H 


4 


C  + 


A^ 


Testing  indicated  that  A-==  2X,  so  (39)  is  simplified  to 

r  . 


A?   =    A 


-t- 


+  »     Cc+0 


which  further  reduces  to 


A^„  =  0.?  A  +■  O.b  A* 


(36) 


(37) 


(33) 


(39) 


(UO) 


(10.) 
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For  simplicity  of  application,  Eq.  (Ip.)  was  further  simplified  to 

£Z0  =  A  +  0.5AZ5-  #  (U2) 

which  was  the  equation  utilized.  However,  testing  indicated  that  Eq. 
(hi)   may  have  given  better  results  than  the  simplified  form  of  Eq.  (U2). 

2.  Testing  of  the  prediction  equation. 

The  prognostic  period  of  12  hours  was  chosen  to  test  the  prediction 
equation.  Four  successive  12-hr  intervals  were  prognosticated  starting 
with  the  0030Z  chart  of  11;  November,  195$ •   Since  the  purpose  of  this 
paper  is  to  attempt  improvement  in  graphical  prognosis  by  the  inclusion 
of  terrain  effects,  the  prediction  equation  was  aoplied  in  two  different 
forms.  First,  a  prognosis  was  made  applying  Eq.  (3U)j  and  a  second  prog- 
nosis using 

^|Vz0+CK+J(J  =  -|-j[ze+ch+J^),?0-Zo+chtJ<Jij3) 

which  is  simply  (3l±)  neglecting  vertical  motion  induced  by  terrain.  With 
minor  differences,  Eq.  (U3)  is  the  prediction  equation  derived  by  Reed  \hj. 
In  applying  (3U)>  the  following  steps  were  carried  out: 

(a)  The  current  "ZQ   chart  was  constructed  by  conventional  graphical 
methods,  using  a  grid  distance  of  600Km. 

(b)  The  current  thickness  chart  was  multiplied  by  the  constant  C 
and  added  to  the  1Q   chart.  In  practice,  using  a  value  of  0.75  for  C  in 
this  step  was  found  to  give  better  results  than  the  derived  value  of  0.9, 
and  consequently  this  value  was  used. 

(c)  Considering  k  =  1,  the  coefficient  of  H  in  the  quantity  F  was 
evaluated  to  be  7.35  x  10"^.  A  terrain  contour  chart  was  multiplied  by 
this  factor  to  produce  the  quantity  F,  and  isolines  of  F  were  added  to 
contours  of  J(cp).  The  chart  of  F  4-  J  (if)  (Fig.  2)  was  added  to  contours 
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Fig.  2 
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produced  in  step  (b)  to  produce  the  advective  field  ( Z0  -+-  Ch  -f  J(cf)  ■+-  F). 

(d)  The  current  ZQ   chart  was  graphically  subtracted  from  the  ad- 
vective field  to  produce  the  conserved  quantity  (ZQ   -  ZQ  +  Ch  +  J(<f)  -f  F) . 

(e)  Charts  (c)  and  (d)  were  superimposed,  and  isolines  of  the  con- 
served quantity  advected  geostrophically  with  the  advective  field  for  the 
prognostic  period  of  12  hours.  A  difficulty  arose  in  the  advection  process 
along  the  west  coast  of  the  United  States  and  Canada  due  to  the  extreme 
gradients  of  the  advective  field  caused  by  the  high  coastal  mountain 
ranges.  It  was  reasoned  that  where  the  terrain  gradient  was  very  large 
the  horizontal  wind  would  not  be  fully  translated  into  upslope  motion, 
and  therefore  result  in  lesser  vertical  velocities.  Thus  the  scalar  multi- 
plier k  in  Eq.  (32)  may  be  appreciably  less  than  unity  in  regions  of  such 
strong  terrain  gradients.  The  amount  of  reduction  of  k  was  determined  by 
cases  tested.  This  reduction  factor  may  be  taken  into  account  by  appro- 
priately reducing  the  gradients  of  F  where  the  gradient  is  very  large. 

The  final  form  of  F  is  given  in  Fig.  2. 

(f)  The  advected  isolines  of  (ZQ  -  ZQ-\r   Ch  +  J(q>)  +  F)  were  sub- 
tracted graphically  from  the  initial  isolines  of  the  same  quantity  to 
produce  isolines  of  A. 

(g)  Half  of  the  12-hr  500-mb  height  change  was  added  graphically 
to  A  to  produce  the  AZQ  chart.  In  practice,  any  suitable  method  of  prog- 
nosticating the  500- mb  chart  could  be  employed,  but  for  this  test  the 
actual  500-mb  height  change  was  used  in  order  that  the  prediction  equa- 
tions, Eqs.  (3U)  and  (k3)  f   might  be  more  objectively  compared. 

In  applying  Eq.  (U3),  the  steps  were  identical  to  the  above,  with 
the  exception  that  in  step  (c)  only  the  contours  of  JT  (cp)  were  added. 

The  results  of  the  above  tests  are  displayed  in  Table  1,  which  com- 
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TABLE  1 


Correlation  coefficient 

:  A 

0.832 

0.8U7 

0.710 

0.790 

B 

0.707 

0.673 

0.580 

0.630 

C 

0.6U 

0.851 

0.790 

0.812 

D 

0.866 

0.870 

0.92U 

0.8U8 

Root-mean-square  error 

:  A 

9iift 

116ft 

103ft 

99ft 

B 

115ft 

172ft 

169ft 

160ft 

G 

85ft 

102ft 

88ft 

89ft 

D 

67ft 

7Uft 

53ft 

69ft 

A  :  Prognostication  including  terrain  factor;  computations  for 
area  west  of  96W,  i.e.,  Rocky  Mountain  area. 

B  :  Prognostication  excluding  terrain  factor;  computations  for 
area  west  of  96W. 

C  :  Prognostication  including  terrain  factor;  computations  for 
entire  area  covered  by  grid. 

D  :  Prognostication  including  terrain  factor;  computations  for 
area  east  of  96W. 

1  :  12 -hr  prognostic  period  commencing  0030Z  lU  November,  1955. 

2  :   "      "  "  "     1230Z  "     "  "  . 

3  :   "      "  "  "     0030Z  15    "  "  . 
h   :   "      "  "  "     1230Z  "     "  "  . 
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pares  height  changes  as  prognosticated  by  the  two  methods  to  the  observed 
height  changes.  The  statistical  parameters  were  computed  using  a  I4.8  point 
square  grid,  which  enclosed  the  area  bounded  by  points  28N,  I2I4.W5  55N, 
I36W5  53N,  65W)  and  27N,  79W.  Figs.  3,  h,   and  5  illustrate  height  changes 
for  the  period  1230Z  15  November  to  0030Z  16  November,  1955.  Fig.  3  shows 
observed  height  changes,  Fig.  1;  the  height  changes  as  prognosticated  by 
Eq.  (3k) >   and  Fig.  5  the  height  changes  as  prognosticated  by  Eq.  (Ii3). 
The  grid  encompasses  an  area  slightly  larger  than  the  above  illustrations. 

3.  Conclusions. 

The  inclusion  of  the  terrain  factor  in  this  model  results  in  materi- 
ally improved  prognoses. 

However,  the  process  of  making  a  prognosis  by  this  method  is  rather 
lengthy,  and  whereas  the  inclusion  of  the  terrain  factor  into  the  pre- 
diction equation  adds  no  more  steps  to  the  routine,  (F  may  be  combined 
with  J((p)  on  a  permanent  chart  and  used  for  each  prognosis),  the  advec- 
tive  procedure  is  much  more  delicate.  The  advective  field  has  areas  of 
very  strong  gradient  with  the  isolines  of  the  conserved  quantity  closely 
packed  and  almost  parallel  to  the  advective  field.  Consequently  consider- 
able care  must  be  taken  to'  insure  accurate  advection. 

In  constructing  the  contours  of  F,  terrain  irregularities  were 
smoothed  to  some  extent.  Some  smoothing  is  mandatory,  as  a  minutely 
accurate  F  chart  is  not  reasonable  to  use  when  considering  large  scale 
flow.  Also,  it  is  evident  that  the  smoother  the  F  chart,  the  easier  the 
graohical  constructions  will  be.  Further  testing  might  indicate  an  op- 
timum balance  between  detail  of  the  terrain  features  and  representation 
for  large  scale  flow. 

As  noted  in  Sect.  2,  a  value  of  0.75  was  used  as  a  multiplier  of  h 
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Fig.  k 
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Fig.  5 
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in  producing  the  advective  field  rather  than  the  derived,  value  of  0.9. 
More  extensive  testing  would  doubtless  indicate  the  optimum  value  of 
this  multiplier.  Reed  [k\    found  that  a  value  of  0.5  gave  good  results 
in  the  three  cases  tested  with  his  model. 
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